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Abstract 

New high precision measurements of the Collins and Sivers asymmetries of charged hadrons 
produced in deep-inelastic scattering of muons on a transversely polarised 6 LiD target are 
presented. The data were taken in 2003 and 2004 with the COMPASS spectrometer using 
the muon beam of the CERN SPS at 160 GeV/c. Both the Collins and Sivers asymmetries 
turn out to be compatible with zero, within the present statistical errors, which are more 
than a factor of 2 smaller than those of the published COMPASS results from the 2002 data. 
The final results from the 2002, 2003 and 2004 runs are compared with naive expectations 
and with existing model calculations. 
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1 An introduction to transverse spin physics 

1.1 Historical introduction 

The importance of transverse spin effects at high energy in hadronic physics was first 
suggested by the discovery in 1976 that A hyperons produced in pN interactions exhibited an 
anomalously large transverse polarisation [1]. This effect could not be easily explained. For a 
long time it was believed to be forbidden at leading twist in QCD [2], and very little theoretical 
work was devoted to this field for more than a decade. Nevertheless some important theoretical 
progress for the understanding of single spin asymmetry (SSA) phenomena was done at that 
time [3]. 

This situation changed in the nineties. After the first hints of large single transverse spin 
asymmetries in inclusive 7r° production in polarised pp scattering at CERN [4] and in IHEP [5], 
remarkably large asymmetries were found at Fermilab both for neutral and charged pions [6] . 
The discovery of the EMC collaboration at CERN in the late eighties that the quark spin con- 
tributes only a small fraction to the proton spin [7] caused a renewed interest in the origin of the 
nucleon spin and proposals for new and versatile experiments. In parallel, intense theoretical 
activity was taking place: the significance of the quark transversity distribution, already intro- 
duced in 1979 [8] to describe a quark in a transversely polarised nucleon, was reappraised [9] in 
1990, and its measurability via the Drell-Yan process established. In 1991 a general scheme of 
all leading twist and higher-twist parton distribution functions was worked out [10], and in 1993 
a way to measure transversity in lepton nucleon polarised deep-inelastic scattering (DIS) was 
suggested [11]. On the experimental side, the RHIC-Spin Collaboration [12] and the HELP Col- 
laboration [13] put forward the first proposals to measure transversity. Today transversity is an 
important part of the scientific programme of the HERMES experiment at DESY, of the RHIC 
experiments at BNL, and of the COMPASS experiment at CERN, all presently taking data. An 
experiment to measure transversity is being prepared at JLAB. First results on a transversely 
polarised proton target have been published recently by the HERMES Collaboration [14] and, 
on a transversely polarised deuteron target by COMPASS [15]. 

The COMPASS published results refer to the data taken in 2002. Further data taking 
occurred in 2003 and 2004, and in this paper results from these new data are presented, as well 
as a reanalysis of the 2002 data and the final results for the 2002-2004 data. 

1.2 The Collins mechanism 

To fully specify the quark structure of the nucleon at the twist-two level, the transverse 
spin distributions A T q(x) must be added to the momentum distributions q(x) and the helicity 
distributions Aq(x) [10], where x is the Bjorken variable. For a discussion on notation, see 
Ref. [16]. If the quarks are collinear with the parent nucleon (no intrinsic quark transverse 
momentum k T ), or after integration over k T , these three distributions exhaust the information 
on the internal dynamics of the nucleon. More distributions are allowed admitting a finite k T , 
as we will see in Sections fT31 and fL4l or at higher twist [16-19]. 

The distributions At<? are difficult to measure, since they are chirally odd and therefore 
absent in inclusive DIS. They may instead be extracted from measurements of the single-spin 
asymmetries in cross-sections for semi-inclusive DIS (SIDIS) of leptons on transversely po- 
larised nucleons, in which a hadron is also detected in the final state. In these processes the 
measurable asymmetry is due to the combined effect of A T q and another chirally-odd function 
which describes the spin-dependent part of the hadronisation of a transversely polarised quark 
q into a hadron h. 

The existence of an azimuthal asymmetry in transversely polarised leptoproduction of 
spinless hadrons at leading twist, which depends on a T-odd fragmentation function and arises 
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Figure 1: Definition of the Collins and Sivers angles. The vectors pj,, s and £? are the hadron 
transverse momentum and spin of the initial and struck quarks respectively. 



from final-state interaction effects, was predicted by Collins [11] and is now generally known 
as the Collins effect. It is responsible for a left-right asymmetry in the fragmentation of trans- 
versely polarised quarks. In some models [20], the Collins asymmetry is expected to be largest 
for the leading hadron in the current jet, i.e. the hadron with the highest momentum. 

Assuming the detected hadron to be spinless, and a collinear quark distribution in the 
nucleon, at leading twist the SIDIS cross-section can be written as 

da 2a 2 s sr^ <> f 1 



dx dy dz d 2 p^ dcf) $ Q 4 



E e * ■ + ( x - vf\ ■ x ■ ■ d >^t) + 

q ^ 



+ (1 - y) \S ± \ sin(<^ - • x ■ A T q(x) ■ A°L>J(z,p£)} (1) 

where Bjorken x is Q 2 /[2M(Ei — Ey)], y is the fractional energy of the virtual photon, Q 2 the 
photon virtuality, M the target nucleon mass, and z = E h / (Ei — Ey) the fraction of available 
energy carried by the hadron. The energies E h , E h and Ey are the energies of the hadron, the in- 
coming lepton, and the scattered lepton respectively in the target rest frame system. The hadron 
transverse momentum p£ is evaluated with respect to the virtual photon direction. Referring 
to Fig. [H 4>h and <f> s i are the azimuthal angles of the hadron and of the struck quark spin in a 
coordinate system in which the z-axis is the virtual photon direction, and the x-z plane is the 
lepton scattering plane with positive x-direction along the scattered lepton transverse momen- 
tum. S± is the target spin normal to the virtual photon direction and <f> s is its azimuthal angle 
with respect to the lepton scattering plane. 

The Pt -dependent fragmentation function can be obtained by investigating the fragmen- 
tation of a polarised quark q into a hadron h, and is expected to be of the form 

D T h g {z,p£) = D h q (z,p£) + A° T D h q (z,p£) ■ sm(<f> h - (2) 



where A^D q (z,p^) is the T-odd part of the fragmentation function, responsible for the left- 
right asymmetry in the fragmentation of the transversely polarised quark. 

The "Collins angle" $c was originally defined in [1 1] as the angle between the transverse 
momentum of the outgoing hadron and the transverse spin vector of the fragmenting quark, i.e., 

(p£ x q)-s' 

sin$ c = l : T q) (3) 
\Pj, x q\ \s'\ 

or 

$c = <Ph - 4> s > ■ (4) 
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Since, as dictated by QED, the directions of the final and initial quark spins are related to each 
other by <f> s i = n — S , equation © becomes 



$c = <Ph + 05 - vr . (5) 

By comparing the cross-sections on oppositely polarised target nucleons one obtains from 
expression <HJ the transverse single spin asymmetry 



A h da(S_L) - da(-S_L) . , . 

where the "Collins asymmetry" is 



t = ; j5 : ; ; s : = P±l " ^aw • ^co« ■ sm y> c (6) 

da(S ± ) + da(-S ± ) 



and 



1 - y + y 2 /2 



is the transverse spin transfer coefficient from the initial to the struck quark [16]. 

1.3 The Sivers mechanism 

An entirely different mechanism was suggested by Sivers [21] as a possible cause of the 
transverse spin effects observed in pp scattering. This mechanism could also be responsible for 
a spin asymmetry in the cross-section of SIDIS of leptons on transversely polarised nucleons. 
Sivers conjecture was the possible existence of a correlation between the transverse momentum 
k T of an unpolarised quark in a transversely polarised nucleon and the nucleon polarisation 
vector, i.e. that the quark distribution q(x) in expression CQ) could be written as 

q T (x, k T ) = q(x, k T ) + \S±\ • A q(x, k T ) ■ sin $ 5 (9) 
where the "Sivers angle" 

<$>s = <P q - <Ps (10) 

is the relative azimuthal angle between the quark transverse momentum k T and the target spin 
S ± . 

Under the assumption that the hadron produced in the fragmentation and the fragmenting 
quark are collinear, i.e. that all the hadron transverse momentum originates from the intrinsic 
transverse momentum of the quark in the nucleon (p^ = zhr), the Sivers angle, shown in Fig.[H 
becomes 

$S = <t>h ~<Ps (11) 

and the SIDIS cross-section in leading order QCD is given by 

dxdydzd 2 p$d(f>s Q 4 Y 2 L J 

q{x,p£/z) + \S ± \ sin(® s )ATq(x,p£/z)] ■ D h q {z) . (12) 
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Comparing the cross-sections on oppositely polarised target nucleons, the transverse spin 
asymmetry of expression © becomes 



dcr(S ± ) + d(j(-S±) 



(13) 



where the "Sivers asymmetry" 



, Z q e 2 q -AZ q (x,p£/z)-D h q (z) 

could be revealed as a sin $5 modulation in the number of produced hadrons. 

The existence of the Sivers function requires final/initial state interaction, and an inter- 
ference between different helicity Fock states. In the absence of interactions the Sivers function 
would vanish by time-reversal invariance of QCD (see e.g. Ref. [19]) and indeed it was believed 
for several years that the Sivers function is zero. Recently it was shown however [22-24] that 
these interactions are represented naturally by the gauge link that is required for a gauge invari- 
ant definition of a transverse momentum dependent (TMD) parton distribution, thus the Sivers 
function has become a very important piece in the most fundamental issues of QCD. 

1.4 More general formalism 

A combined treatment of the the Collins and Sivers mechanisms requires a more general 
approach and the introduction of TMD distributions and fragmentation functions [18, 19, 25]. 
Convolution integrals appear in the SIDIS cross-section, whose expression becomes 



da ~ ^ej-|i[l + (l-y) 2 ] x [ g ®Dj+|SJsin($ s )Afo®££ 



+ 



\S ± \ S m(<S> c )xA T q®A°,D h q } . (15) 

In the general expression of the SIDIS cross-section at leading order QCD other terms related 
to different single and double spin azimuthal asymmetries appear. Here they are neglected since 
they are beyond the scope of this paper. The symbol ®, which replaces the products in Eq. © 
and (fT2"l) . indicates the convolution integral 

DF ® FF = J d 2 k T DF(x, k T ) ■ FF(z,p£ - zk T ) (16) 

where DF and FF are generic TMD distributions and fragmentation functions respectively. 
The transverse spin asymmetry is then given by 

A h = da(S ± ) - do-(-Sj_) 

T ~ d(7(S ± ) + dflr(-Si) 

= I S±_ I • D NN Acoii ■ sin + | S± | • A Siv ■ sin $ s (17) 

where the Collins and the Sivers asymmetries are still given by Eq. © and (TbH) when replacing 
the products of the distribution and fragmentation functions with the corresponding convolu- 
tions. In the usual assumption of Gaussian distributions for the parton and the hadron transverse 
momenta in the DF and in the FF, the only effect of the convolution integral in Eq. and 
(fT4l is the presence of a factor which depends on < fc T 2 > and < p£ 2 > (see f.i. [26, 27]). 

Since the Collins and Sivers terms in the transverse spin asymmetry depend on the two 
independent angles $c and $5, measuring SIDIS on a transversely polarised target allows the 
Collins and the Sivers effects to be disentangled and the two asymmetries can separately be 
extracted from the data. 
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2 The COMPASS experiment 

2.1 Physics objectives 

The COMPASS experiment was proposed to CERN in 1996 to investigate hadron struc- 
ture and hadron spectroscopy by carrying on a number of key measurements using both hadron 
(7r, K and protons) and muon high-energy beams. Apart from a short pilot run in 2004 with 
190 GeV/c pions to measure the pion polarisability via Primakoff scattering on high-Z targets, 
the experiment in so far has focused on the investigation of the spin structure of the nucleon 
using a 160 GeV/c ji + beam and a polarised deuteron target. In the longitudinal target spin 
mode, the main goal of the experiment is the measurement of AG/G [28], the polarisation of 
the gluons in a longitudinally polarised nucleon, but very precise A\ data are also collected [29] . 
In about 20% of the running time, the target polarisation was set along the vertical direction, 
orthogonal to the beam axis, and transverse spin effects were measured, which are the subject 
of this paper. 

2.2 The experimental set-up 

The COMPASS spectrometer [30] has been set up in the CERN SPS North Area. 

The layout of the spectrometer which was on the floor in 2003 is shown in Fig. El To 
combine large geometrical acceptance and broad dynamical range, the spectrometer comprises 
two magnetic stages, and is made up of a variety of tracking detectors, a fast RICH, two hadronic 
calorimeters, and provides muon identification via filtering through thick absorbers. 

The first stage is centred around the spectrometer magnet SMI (1 Tm bending power) and 
the design acceptance is about ±200 mrad in both planes, to fully contain the hadrons of the 
current jet. The second stage uses the spectrometer magnet SM2 (operated at a bending power 
of 4.4 Tm), located 18 m downstream from the target, and the acceptance is ±50 and ±25 mrad 
in the horizontal and vertical plane respectively. 

The design of detector components, electronics and data acquisition system allows to 
handle beam rates up to 10 8 muons/s and about 5-10 7 hadrons/s. The triggering system and the 
tracking system of COMPASS have been designed to stand the associated rate of secondaries, 
and use state-of-the-art detectors. Also, fast front-end electronics, multi-buffering, and a large 

Hadron 



Hadron absorber 2 




10 20 30 40 50 [m] 



Figure 2: Top view of the lay-out of the spectrometer for the COMPASS experiment in 2003. 
The labels and the arrows refer to the major components of the trigger and PID systems. The 
thin vertical lines represent the tracking detectors. 
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and fast storage of events are essential. 

The experiment has been run at a muon energy of 160 GeV. The muons originate from the 
decay of n and K mesons produced by the 400 GeV proton beam on a primary beryllium target. 
The beam is naturally polarised by the weak decay mechanism, and the beam polarisation 
is about -80%. The beam polarisation contributes to the transverse spin dependent part of the 
cross-section only by higher twist effects, which are not considered in the leading-order analysis 
of this paper. The fi + intensity is 2 ■ 10 8 per spill of 4.8 s with a cycle time of 16.8 s. The beam 
profile presents a Gaussian core and a large non-Gaussian tail due to halo muons. The beam has 
a nominal energy of 160 GeV and is focused at the target centre, with a spread of 7 mm (r.m.s.) 
and a momentum spread of a p /p = 0.05 for the Gaussian core. The momentum of each muon is 
measured upstream of the experimental area in a beam momentum station consisting of several 
planes of scintillator strips with 3 dipole magnets (30 Tm total bending power) in between. The 
precision of the momentum determination is typically Ap/p < 0.003. 

Due to the major problems and delay in the construction of the new large-acceptance 
COMPASS polarised target (PT) magnet, the experiment has utilised in so far the magnet from 
the SMC experiment, which has a similar design and identical magnetic properties, but a smaller 
bore (26.5 cm diameter). The resulting angular acceptance is reduced, going from ±170 mrad 
at the downstream end to ±69 mrad at the upstream end of the target. The polarised target sys- 
tem [31] consists of two oppositely polarised target cells (upstream u, downstream d), 60 cm 
long each and 3 cm diameter, so that data are collected simultaneously for the two target spin 
orientations. The PT magnet can provide both a solenoidal field (2.5 T) and a dipole field for 
adiabatic spin rotation (up to 0.5 T) and for the transversity measurements (set at 0.42 T). Cor- 
respondingly, the target polarisation can then be oriented either longitudinally or transversely to 
the beam direction. The target is cooled to temperatures below 100 mK by a 3 He- 4 He dilution 
refrigerator. When operated in the frozen spin mode, temperatures of ~50 mK are reached. The 
polarisation was lost at rate of 0.4 - 1.0%/day in the 0.42 T dipole field and 0.05 - 0.1%/day in 
2.5 T solenoid field. The use of two different target materials, NH 3 as proton target and 6 LiD as 
deuteron target, is foreseen. Polarisations of 90% [32] and 50% [30] have been reached, respec- 
tively. In so far only 6 LiD has been used as target: its favourable dilution factor of ~0.4 is of the 
utmost importance for the measurement of AG. The dynamical nuclear polarisation system can 
polarise the target only at 2.5 T, i.e. in the solenoidal field. To run in the transverse polarisation 
mode, the spins are frozen and rotated adiabatically by first lowering the longitudinal field to 
0.5 T, then rotating the magnetic field to the vertical direction with the help of the dipole coils. 
The polarisation values are obtained from corresponding measurements done at the beginning 
and at the end of each data taking period, with the polarised target field set back at 2.5 T, and 
take into account the relaxation time of the target polarisation. 

To match the expected particle flux in the various locations along the spectrometer, COM- 
PASS uses very different tracking detectors. The small area trackers consist of several stations 
of scintillating fibres, silicon detectors, micromegas chambers [33] and gaseous chambers using 
the GEM-technique [34]. Large area tracking devices are made from gaseous detectors (Drift 
Chambers, Straw tubes [35], and MWPC's) placed around the two spectrometer magnets. 

Muons are identified in large-area Iarocci-like tubes and drift tubes downstream of muon 
absorbers. Hadrons are detected by two large iron-scintillator sampling calorimeters, installed 
in front of the absorbers and shielded to avoid electromagnetic contamination. The charged 
particles identification relies on the RICH technology [36]. In this paper we have not utilised 
the information of the RICH, and give results for non-identified hadron asymmetries only. The 
asymmetries for RICH-identified hadrons (pions and kaons) will be the subject of a separate 
paper. 
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Figure 3: General architecture of the DAQ system. 



The trigger [37] is formed by a combination of signals indicating the presence of a scat- 
tered muon at a given angle or in a given energy range. In most DIS events (Q 2 > 1 (GeV/c) 2 ), 
the scattered muon is identified by coincidence signals in the trigger hodoscopes, that mea- 
sure the projection of the scattering angle in the non-bending plane and check its compatibility 
with the target position Several veto counters installed upstream of the target are used to avoid 
triggers due to halo muons. In addition to this inclusive trigger mode, several semi-inclusive 
triggers select events fulfilling requirements based on the muon energy loss and on the presence 
of a hadron signal in the calorimeters. The acceptance is further extended toward high Q 2 values 
by the addition of a standalone calorimetric trigger in which no condition is set for the scattered 
muon. 

A complete description of the spectrometer can be found in Ref . [30] . 

2.3 Data taking and off-line system 

The data acquisition system [38] was designed to read the 250000 detector channels 
with an event rate of up to 100 kHz and with essentially no dead-time. This required a full 
custom design of the readout electronics (mounted directly on the detectors) and of the readout- 
driver modules (named CATCH and GeSiCA), which perform local event building and trigger 
distribution to the front-end boards (Fig. fSJ). The trigger control system (TCS) performs trigger 
distribution and synchronisation of the time-to-digital converters to better than 50 ps. The data 
collected during the 4.8 s spills is stored in 32 GByte buffers. The event building is performed 
on high performance Linux PCs connected to the read out system via Gigabit Ethernet. The 
data are grouped in "runs" corresponding to about 100 consecutive spills in 2003 and 200 in 
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Figure 4: Schematic view of the off-line system and reconstruction and analysis flow. 

2004. General information on the run is extracted and stored in the meta-data tables of an Oracle 
database. The files are then sent to the Central Data Recording facility in the CERN computer 
centre in parallel multiple streams over a dedicated optical fibre network at an average speed of 
70 MBytes/s. 

Once at the computer centre, the files are registered in the name space of CASTOR (the 
CERN hierarchical storage management system): from this moment onward, CASTOR controls 
thoroughly the events data handling (copy to tape, managing of the disk space), while an Oracle 
RAC database system is in charge of translating high-level requests of data into file requests. 
The huge amount of data of about 350 TBytes/year is reconstructed at the CERN computer 
centre, requiring a computing power of about 200k SPECint2000. The event reconstruction is 
carried on by CORAL, the COMPASS reconstruction and analysis framework, a fully object- 
oriented program with a modular architecture written in C++, which provides interfaces for 
the event reconstruction algorithms and insulation layers to access the data and for external 
pluggable packages [39]. The reconstruction is carried out in parallel, by some 600 jobs running 
on the CERN batch system. CORAL decodes the data, reconstructs tracks and vertexes, and 
performs the particle identification making use of the alignment and calibration data describing 
the apparatus, which are stored as time/version-dependent information in a MySql database. 
The reconstructed data is output in a proprietary format to files called Data Summary Tapes 
(DSTs), and in ROOT [40] format to mini-DSTs, which are selectively filtered out from the 
DSTs during production and turn out to have a size of about 1% of the original RAW data. 
DSTs and miniDSTs are stored also centrally on tape, under CASTOR. The physics analysis 
is performed on the mini-DSTs, replicated in the different institutes, by means of PHAST, the 
COMPASS framework for the final data analysis [39]. Fig. |4] depicts the reconstruction and 
analysis system to which the data flow after they are stored centrally. For the simulations, the 
Monte Carlo program COMGeant [39], based on GEANT 3 and LEPTO 6.5.1, has been written 
and used extensively. 
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In the first three years of data taking, from 2002 to 2004, in which the experiment ben- 
efited of about 12 ■ 10 5 spills, COMPASS collected 30 billion events, corresponding to a total 
data sample of more than 1 PByte. Thanks to continuous work on alignment and improvement 
of the reconstruction code, these data have been processed several times. About 20% of these 
data have been taken in the transverse target spin mode. 

2.4 Principle of the measurement 

As explained in Section fTTTl the Collins and the Sivers asymmetries can be estimated from 
the SIDIS data by measuring, in a given kinematical bin defined by x, z, or p\, the number of 
events N + and N~ collected on oppositely polarised target nucleons, and fitting the measured 
asymmetries (iV + — iV~)/(iV + — N~) either with a sin$ c or a sin$ s dependence. Since 
the spectrometer acceptance is somewhat different for the two target halves, the cross-section 
asymmetries cannot be obtained from a direct comparison of the number of events collected in 
the u — and d— cell. It is necessary to compare the number of events collected in each cell with 
the two target spin orientations, and the two-cell system allows to reduce the systematic effects. 

The orientation of the target polarisation cannot be flipped fast in the COMPASS target 
system. The holding field is 0.42 T, and transverse to the beam, thus it is out of question to 
flip the target spin by inverting the magnetic field because the acceptance would change. Po- 
larisation reversal is done by exchanging the microwave frequencies for w— and d— cells, and 
to achieve ~50% polarisation one needs two days. Therefore, once the target is polarised, data 
are taken for several days (typically 5) before a polarisation reversal is done, and to evaluate 
an asymmetry one compares the number of events collected almost a week apart. For this rea- 
son, the transverse spin measurement usually was carried out at the end of the run, when the 
spectrometer was fully operational and functioning smoothly, and great care was taken not to 
intervene on any detector to avoid changes in efficiency. 

The typical cycle (a data taking period) consisted therefore of 5 days of measurement in 
one configuration (yielding f.i. and N~[), 2 days to change the target polarisation, and 5 
more days of data taking in the opposite configuration (yielding N~ and N£). At this point, 
one can estimate the cross-section asymmetry. A straightforward procedure was used for the 
published 2002 data analysis [15], namely to derive one asymmetry from iV+ and N~, a second 
asymmetry from Nj~ and N d , check their consistency and then average the two. A different 
method, which uses at once the four numbers and is less sensitive to differences in acceptance, 
has been used in this analysis, and is described in Section l4.3.2l 

In this paper final results are presented for all the data collected with the deuteron target in 
the three years 2002-2003-2004. Data were collected during two periods in 2002, one period in 
2003, and during two periods in 2004. The focus is on the data collected in 2003 and 2004, and 
most of the next section will be devoted to the analysis of these data, which is slightly different 
from the analysis of the 2002 data (already published), thanks to various improvements in the 
analysis code. The new chain of analysis has also been used to reanalyse the 2002 data: the 
results are perfectly compatible, as will be shown in Section RTTl Also, the new 2003-2004 data 
are in excellent agreement with the 2002 data. At the end, overall results will be given for the 
whole data set, 2002-2004. 

3 Data analysis I - event reconstruction and selection 
3.1 Event reconstruction 

In the event reconstruction, a track reconstructed before the target is assumed to be an 
incoming muon if it is reconstructed in the scintillating fibres and silicons upstream of the target, 
its momentum is reconstructed in the BMS, and the track time is within 3 standard deviations 
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of BMS and trigger time. If several valid BMS-tracks are compatible with the track time, a 
backtracking algorithm is used to resolve ambiguities. 

The scattered muon candidates are defined as the positively charged outgoing tracks with 
a momentum larger than 1 GeV/c, going through SMI, and their extrapolation at the entrance 
and at the exit of the target is within 5 cm from the target axis. In addition, for all triggers based 
on the hodoscope information, the track must be compatible with the hodoscope hits as given 
in the trigger matrix. In the case of a calorimetric trigger, a minimal number of hits is required 
in the muon walls and the amount of material traversed in the spectrometer must be larger than 
66 and 74 radiation lengths for tracks reconstructed in the first and in the second spectrometer 
stage respectively. 

The muon interaction point (the so-called "primary vertex") is defined by one beam par- 
ticle and one scattered muon. If in the event there are more than one beam particles and/or 
scattered muons, several vertexes may be reconstructed. The distance of closest approach be- 
tween the scattered muon and beam track must be less than 10 a (about 2-3 mm). The position 
of the vertex along the spectrometer axis z vtx ^ is given by the average of the distance of closest 
approach of all tracks in the vertex, while in the orthogonal plane the coordinates x vtx and y vtx 
are defined from the beam track at z vtx . The tracks belonging to the vertex are selected using a 
Kalman fit. The "best primary vertex", used in the following steps of the analysis, is defined as 
the one with the maximum number of tracks and, if the number of tracks is the same, the one 
with smaller vertex x 2 - 

Only the events with at least one primary vertex reconstructed with at least one more 
outgoing track, and Q 2 > 1 (GeV/c) 2 , which are about 1% of the initial raw event sample, are 
written on miniDST and used in the physics analysis. 

The event reconstruction was almost the same for the 2003 and 2004 data, and is very 
similar to that used for the 2002 data. Each year an improvement in the ratio of reconstructed 
events over the useful beam was achieved, due both to additional detectors in the spectrometer 
and to the better tuned reconstruction and analysis software. 

3.2 Data quality checks 

To monitor the performances of the apparatus before and after the target spin reversal (i.e. 
on the two "sub-periods") in each data taking period, the time stability of many distributions 
has been checked dividing the whole sample into runs or clusters of neighbouring runs. 

Using the histograms produced during the event reconstruction, the detector performance 
stabilities were scrutinised looking at the stability of the shape of the hit distributions in the 
about 360 detector planes. The time stability of the detector and reconstruction efficiencies was 
checked looking at: 

- the number of clusters per plane and per event, 

- the mean number of tracks per event, 

- the mean number of track segments in the different spectrometer regions per event, 

- the mean number of primary vertexes per event, 

- the mean number of secondary vertexes per event. 

Using the miniDST events, the stability was checked monitoring run per run 

- the number of reconstructed K° per primary vertex, 

- the reconstructed K° mass distribution, 

- the energy measured in the two hadronic calorimeters HCAL1 and HCAL2, 

^ the used reference system is denned to have a right handed frame with the z axis along the spectrometer axis 
(i.e. the nominal beam direction) the y axis in the vertical direction pointing upwards, and the origin at the 
centre of the u-cell. 
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Figure 5: nn invariant mass in the K° region as a function of time from August 13 to August 
19,2004. 



- the x vtx and y vtx distributions in the two cells, 

- the vertex x 2 distribution. 

Also, the time stability of the distributions of several kinematical observables (like x, Q 2 , y, the 
azimuthal angles and the momenta of the scattered muon and of the hadrons) was investigated 
in detail. 

As an example, the mean nix invariant mass in the K° region from the data collected in 
the first sub-period of the 2004 run is shown in Fig.|5]as a function of time. 

Runs showing some instability were not used for the physics analysis. The runs rejected 
by this criterion are 28 (over a total of 458) for the 2003 data, and 44 (over 462) for the 2004 
data, corresponding to about 5% and 4% of the initial raw event sample. 

3.3 DIS events selection 

To better define the DIS events, more refined cuts were applied. In particular: 
1. primary vertex: if more primary vertexes were reconstructed in one event, the best pri- 
mary vertex was selected. 

The vertex had then to be inside the target. A radial cut r vtx < 1.3 cm was applied on the 
distance of the vertex from the target axis. The event was then accepted only if z vtx was 
inside one of the two target cells (-100 cm < z vtx < -40 cm or -30 cm < z vtx < 30 cm), 
and assigned correspondingly to the u- or to the J-cell. 

The distribution of the primary vertex z-coordinate for the final sample is shown in Fig.|6l 
The increase in the number of events with z vtx is due to the increase in geometrical ac- 
ceptance going from the upstream to the downstream end of the target, as explained in 
Section |2~21 and is very well reproduced by the Monte Carlo simulation. The two target 
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Figure 6: Distribution of the primary vertex z-coordinate for the final sample. The dashed his- 
togram indicates the events selected by the z vtx cut. 

cells are clearly separated. One can also notice the continuum of events produced on the 
helium bath, as well as a sample of events produced in the Al window, at z= 50 cm, 
which seals the vessel of the PT magnet. 

2. incoming muon: the cut Pbeam < 200 GeV/c was applied on the momentum of the incom- 
ing particle. A check was also performed on the quality of the beam track, and a safe cut 
on x 2 /^ was applied. 

To ensure an identical beam intensity in both target cells, i.e. a nearly identical luminos- 
ity, the beam track projection at the entrance and at the exit of the target region (e.g. at 
z= —100 cm and z= 30 cm) had to be inside the target region as defined for the primary 
vertex. The momentum distribution for the final data sample is shown in Fig. [7] 

3. scattered muon (/i'): as well as for the incident muon, a suitable cut was applied on the 
X 2 / v of all reconstructed tracks from the vertex, including the scattered muon candidate. 
Other outgoing tracks were flagged as muon candidate if they had hits in one of the two 
Muon Walls. To achieve a clean muon identification, the amount of material traversed in 
the spectrometer had to be larger than 30 radiation lengths. Only events with one and only 
one muon candidate entered the following steps of the analysis. 

In addition, standard DIS cuts Q 2 > 1 (GeV/c) 2 , mass of the final hadronic state W > 5 GeV/c 2 , 
and 0.1 < y < 0.9 were applied. All these cuts reduced the number of miniDST events by a 
further 45%. 

3.4 Hadron identification 

All the outgoing particles not flagged as // candidates were assumed to be hadrons. 
Only particles with the last measured coordinate after the first spectrometer magnet were 
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Figure 7: Momentum distribution of the reconstructed incoming muons for the final sample of 
events. 

used, to reject tracks reconstructed in the fringe field on SMI which have a poorer momentum 
resolution. 

The following requirements had to be satisfied in order to identify the particle with a 
hadron: 

1. the amount of material traversed in the spectrometer had to be smaller than 10 radiation 
lengths; 

2. the particle must not give a signal in any of the two hadronic calorimeters HCAL1 and 
HCAL2. 

3. if the particle gave a signal in one calorimeter, the measured energy had to be E HC ali > 
5 (2003 data) or 4 (2004 data) GeV, or E HCAL2 > 5 (2003 and 2004 data) GeV. The 
correlation between the energy measured in HCAL and that measured by the spectrometer 
is shown in Fig. [8] 

4. the transverse momentum of the particle with respect to the virtual photon direction had 
to be larger than 0.1 GeV/c. 

The first requirement reduced the muon contamination, the second and the third the muon 
and the electron contamination. The last cut was introduced to assure a good resolution 
in the measured azimuthal angle. 
Two different asymmetries have been measured. The "all hadron" asymmetries were evaluated 
using all the particles identified as hadrons carrying a fraction of available energy z > 0.2. For 
the "leading hadron" analysis, only the events with at least one hadron with z > 0.25 were used, 
and the hadron with largest z was defined as leading. In the case the energy of the reconstructed 
leading hadron was higher than the missing energy evaluated from all the charged reconstructed 
hadrons, no further test was performed. If this was not the case, the event was accepted only 
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Figure 8: Correlation between the energy measured in HCAL1 (left) and HCAL2 (right) and 
the energy measured by the spectrometer for the 2004 data. 

if in the hadron calorimeters no particle with energy larger (within the energy resolution) than 
that of the leading hadron was detected. 

After these cuts, the number of events with a leading hadron amounted to 22% of the 
number of miniDST events. 

4 Data analysis II - kinematical distributions and asymmetry evaluation 
4.1 Final samples of events 

In table[H the final statistics used for the asymmetry evaluation is given for all the periods 
both for the positively and negatively charged "leading hadron" and the "all hadron" sample. 



Table 1 : Final statistics used for the asymmetry evaluation 



Year 


Period 


Leading hadron sample 


All hadron sample 






positive hadrons 


negative hadrons 


positive hadrons 


negative hadrons 


2002 


1 


0.48 • 10 6 


0.38 • 10 6 


0.71 • 10 6 


0.59 • 10 6 


2002 


2 


0.32 • 10 6 


0.26 • 10 6 


0.48 • 10 6 


0.40 ■ 10 6 


2003 


1 


1.68 -10 6 


1.33 • 10 6 


2.46 • 10 6 


2.03 • 10 6 


2004 


1 


1.44- 10 6 


1.13 - 10 6 


2.12- 10 6 


1.74- 10 6 


2004 


2 


1.87 • 10 6 


1.47 • 10 6 


2.75 • 10 6 


2.26 ■ 10 6 



4.2 Kinematical distributions 

The distributions of some kinematical quantities from the final 2004 leading hadron 
events are shown in Figs. |5]to[TTJ The plots have been produced after all the cuts described 
in the previous section, if not specified differently. 

The Q 2 vs x scatter-plot and its projections are shown in Fig.|9j 

The x-z correlation (without the z cut), the x-y correlation, and the z-p\ correlation 
(without the z and the p\ cuts) are shown in Fig. |T0l while the y, p\ (without the p\ cut) and 
the z (without the z cut) distributions are shown in Fig.fTTI 
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Figure 9: Scatter-plot of Q 2 vs x, and the corresponding Q 2 - and x-distributions for the 2004 
positive plus negative leading hadron sample. 

4.3 Asymmetry evaluation 

4.3.1 Binning 

The Collins and Sivers asymmetries were evaluated as function of x, p?, and z dividing 
the corresponding kinematical range in bins (with variable width, in order to have a comparable 
statistics in each of them), and integrating over the other 2 variables. No attempt has been 
done to extract the asymmetries in a 2- or 3-dimensional grid. In total, the asymmetries where 
evaluated in the 9 x-bins, 9 p^-bins, and 8 ,2 -bins for the all hadron samples (7 for the leading 
hadron, because of the z > 0.25 cut) given in Tabled 

Tables with the mean values of Q 2 , z, p\ and y in all the x bins are available on HEP- 
DATA [41]. As an example, Fig. fT2l shows the mean values of Q 2 (left), z, p\ and y (right) in 
the different x bins for the all positive hadron sample from 2003 and 2004 data. 

4.3.2 Raw asymmetry evaluation 

The Collins and Sivers asymmetries have been evaluated separately, in each kinematic 
bin, for each data taking period, and for positive and negative hadrons. 

Using Eq. ([TJ and (|T2b . the $ distribution of the number of events for each cell and for 
each polarisation state can be written as 

= *uW ■ (1 ± 4k *n*i)> (18) 



17 



0.8 



0.6 



0.4 



0.2 



icr- 




* fx 2 
Oh 



0.2 0.4 0.6 0.8 



Figure 10: Scatter-plots of z vs x, y vs x, and z vs p 1 ^ for the 2004 all positive plus negative 
hadron samples. 



where j = C, S, and F is the muon flux, n the number of target particles, a the spin averaged 
cross-section, and aj the product of angular acceptance and efficiency of the spectrometer. The 
index k = u,d refers to the target cell. The ± signs indicate the two target spin orientations: 
+ means spin up, — spin down. Here the Collins and Sivers angles &j are evaluated using the 
expressions given in Sections fL2l and n~3l always assuming target spin up, at variance with what 
was done in [15]. 



0.20 < z < 0.25 
0.25 < z < 0.30 
0.30 < z < 0.35 
0.35 < z < 0.40 
0.40 < z < 0.50 
0.50 < z < 0.65 
0.65 < z < 0.80 
0.80 < z < 1.00 



0.10 <v h T 
0.20 < v h T 
0.30 < v h T 
0.40 < Vt 
0.50 < Vt 
0.60 < v h T 
0.75 < v h T 
0.90 < v% 
1.30 <j4 



0.003 < 
0.008 < 
0.013 < 
0.020 < 
0.032 < 
0.050 < 
0.080 < 
0.130 < 
0.210 < 



x < 0.008 
x < 0.013 
x < 0.020 
x < 0.032 
x < 0.050 
x < 0.080 
x < 0.130 
x < 0.210 
x < 1.000 



< 0.20 GeV/c 

< 0.30 GeV/c 

< 0.40 GeV/c 

< 0.50 GeV/c 

< 0.60 GeV/c 

< 0.75 GeV/c 

< 0.90 GeV/c 

< 1.30 GeV/c 



Table 2: Final binning in x, z, and p^. 
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Figure 11: Upper plots: distributions of y and (without the cut) for the 2004 all positive 
plus negative hadron sample. Lower plots: ^-distribution (without the z cut) for the same events 
(left) and for the 2004 leading positive plus negative hadron sample (right). 

The asymmetries e^ k are 

e c,k = c c,k ' A Coii , e S k = Cg k ■ A Siv , (19) 



where 



4,k = f-\PT, k \-D N N, 4j, = f-\P±\. (20) 



NN 



The factor / is the polarised target dilution factor, P^ k the deuteron polarisation, and D 
the transverse spin transfer coefficient from the initial to the struck quark given in Eq. ©. The 
absolute values of the target polarisation Pt in the two cells and for the two spin orientations in 
the same cells differed up to about 7%. 

In order to estimate e^ k from the measured number of events a method, which in the 
following will be called "ratio product method" (RPM), has been used which minimises the 
effect on the estimated asymmetries of possible acceptance variations. 

In each data taking period, the "ratio product" quantities [15] 

)( )) "ivr„ft) N- d ^r ] ~ c ' s <21) 
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Figure 12: Mean values of Q 2 (left), z, p% and y (right) in the different x bins for all positive 
hadrons, 2003 and 2004 data. 

were evaluated from the number of events in the two cells. For small values of the quantities 

4,k k is: 



;i + s>- sin$ i)-( 1 + sV sin< ^ 



where 



1 - e i)U ■ sin $ d ) ■ (1 - e jd ■ sin $,) 
C F -C aJ -(l + Af -sin^) (22) 



A? = 4u + 4a + + e", = 4- < e,- > , (23) 



■ F d Hu(®i) ■ a i,d(®i 



The raw asymmetries A™ and A™ were evaluated by fitting respectively the quantities A c (&c) 
and A s (&s) with the functions p ■ (1 + A™ ■ sin & c ) and p • (1 + A™ ■ sin $5), where p is a 
free parameter. 

In this method the only requirement to avoid systematic errors due to acceptance effects is 
that the factor C a j does not depend on $ 3 . This is surely true, in particular, under the reasonable 
assumption that the ratio Oj U (Q>j)/ aj d ($j) before the polarisation reversal be equal to the 
corresponding ratio aj u (&j)/ a+ d ($j) after the reversal in each $j bin. In this case, C a j is 
equal to 1 for all values. Also, with the cuts applied on the incoming beam it is expected that 
Cp — 1 and indeed the fits to A™ give p values always compatible with 1. 

In addition to the fact that this estimator has "soft" requirements on the acceptance sta- 
bility, it is independent of the relative luminosity and combines all the data from the two target 
cells. A further advantage of the use of the RPM, is that at first order (for small values of the 
involved asymmetries) all spin-independent effect, e.g. Cahn asymmetry, are factored out. 

Concerning the $j binning, the interval (0, 2ir) was divided in 16 bins. Monte Carlo 
studies indicated that the angular resolution (rms) is 0.07 rad, much smaller than the bin size. 

4.3.3 Collins and Sivers asymmetries evaluation 

To extract the Collins and Sivers asymmetries, the target polarisation values given in 
table |3] have been used. 
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Year 


Period/S ubperiod 


Polarisation 






upstream 


downstream 


2003 


1 / 1 


-49.7 


52.8 




2/2 (1 st part) 


49.4 


-42.6 




2 / 2 (2 nd part) 


51.3 


-44.6 


2004 


1 / 1 


50.70 


-43.52 




1 / 2 (1 st part) 


-44.8 


46.0 




1 / 2 (2 nd part) 


-38.6 


40.4 




21 1 (1 st part) 


-46.1 


47.4 




2 / 1 (2 nd part) 


-46.4 


47.4 




2/2 


49.9 


-42.8 



Table 3: Target polarisation values for the 2003 and 2004 data taking. 

The dilution factor / has been taken constant and equal to 0.38 for all data taking periods. 

Using the event kinematics, the mean values of the quantities c^ k = f ■ \P^ k \ ■ Dnn and 
c 5fc = / ' l-^Tfcl nave b een evaluated for each bin and each period of data taking. These values 
have than been used to calculate the Collins and Sivers asymmetries from the corresponding 
raw asymmetries A™. 

5 Data analysis III - Systematic studies 

Given the high statistics of our measurements, a number of systematic studies have been 
performed in order to determine the size of possible systematic errors. 

Extensive tests both to measure false asymmetries and to investigate the stability of the 
physics results were done, for each measured asymmetry, and in each data taking period. They 
are briefly described in the following. 



5.1 False asymmetry studies 

Two different kinds of asymmetries expected to be zero were measured using the final 
event samples, for the Collins and Sivers angles, positive and negative hadrons, leading and all 
hadrons, and in each x, z, "p\ bin. 

The first one was built by splitting the two target cells in two parts and by combining the 
data from the same cell. The mean values of all the resulting asymmetries were found to be 
compatible with zero, as it should. The distributions of the pulls of these asymmetries with re- 
spect to zero, i.e. the values of the resulting asymmetries divided by the corresponding standard 
deviation, were well compatible with Standard Normal distributions, giving no indication of the 
presence of systematic effects. 

The second asymmetry was measured after scrambling the data collected in each period 
and measuring the asymmetry in the standard way. The runs of the two sub-periods before and 
after the polarisation reversal were mixed by labelling as "spin up" every 1st, 3rd, 5th, ... run of 
the two sub-periods and as "spin down" all the others. This selection was truncated in a way, 
that each fake period had the same amount of runs from each real sub-period to ensure roughly 
the same amount of events. As in the previous case, the asymmetries are expected to be zero, 
and the measured asymmetries were compatible with zero. The refined statistical analysis of 
pulls with respect to zero did not give any hint for the presence of systematic effects. 
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5.2 Stability of physics asymmetries vs acceptance and time 

Possible acceptance effects on the physics results from each data taking period were tested 
by checking the compatibility of the asymmetries evaluated after splitting the data according 
to: the location inside the target; the region of the spectrometer in which the scattered muon 
was measured; the trigger of each event. This work was done for both the Collins and Sivers 
asymmetries of positive and negative hadrons, leading and all, and in each x, z, p\ bin. 

The compatibility test consisted in comparing the distribution of the pulls 



with the normal standard distribution. Here A™ are the raw asymmetries evaluated from the 
whole data sample of one specific period k. They are evaluated for positive and negative hadrons, 
Collins and Sivers asymmetries in each x, p\ and z bin, for a total of 2 x 2 x (9 + 9 + 7) = 100 
values in each period. The corresponding variance is a\m . The asymmetries A™ k are the cor- 
responding ones evaluated for all the sub-samples % in which the original data set was divided. 
The variances of A^i are a\ m , and the use of the difference with a 2 Am takes into account almost 
completely all the correlations between A™ k and A™. The compatibility with the normal stan- 
dard distribution was checked also for some specific group of asymmetries, like positive and 
negative hadrons, and Collins and Sivers asymmetries. 

The effect of the different acceptances for events with the primary vertex in the two target 
cells was investigated, evaluating the physics asymmetries separately for the events with the 
primary vertex in the inner half and in the outer half of the two target cells. Combining events 
from only the outer halves (upstream part of cell u and downstream part of cell d) as well as 
from the inner halves only (downstream part of cell u and upstream part of cell d) gave results 
compatible to each other and to the ones using the full sample. 

The asymmetries were also evaluated dividing the data samples accordingly to the az- 
imuthal angle "Jy of the scattered muon in the laboratory system, namely in the regions < 
<Ey < 7r/2, 7r/2 < 3y < ir, n < <Jy < 3n/2, and 3n/2 < 5y < 2p. 

As in the previous case, the pull distributions turned out to be centred around zero, and 
in very good agreement with the standard distribution. The RMS were statistically compatible 
with 1, as expected in the absence of systematic effects. 

The physics asymmetries have been evaluated splitting the data in different samples ac- 
cording to the different triggers. The results turned out to be statistically compatible, giving 
once more no evidence for systematic effect. 

A further test was done to investigate the stability in time of the physics asymmetries. 
To do this, each of the two sub-periods entering the standard asymmetry calculation was split 
in these 9 groups of subsequent runs (each corresponding to 10 hours of data taking). The 
asymmetries were then calculated using the full data set of one sub-period and all the 9 groups 
of runs of the other sub-period separately. In total 54 asymmetries were evaluated from the 
2003 and 2004 data. This method is known to be quite powerful to single out significant time 
dependencies within single periods. As for the previous test, the pulls relative to the mean 
asymmetry for each time slot region were calculated, but no evidence for systematic errors was 
found in their distributions. 

A further test was done to check the possible effect of z and p\ acceptance. The asym- 
metries were evaluated also in 4 complementary z and p\ bins: (z < 0.5, < 0.5), (z < 0.5, 
p\ > 0.5), (z > 0.5, p\ < 0.5), (z > 0.5, p\ > 0.5) and their weighted means were compared 
with the standard results. The differences between the two sets of measurements are essen- 




(25) 
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tially invisible, much much smaller than the statistical errors, again as expected in absence of 
systematic effects. 



5.3 Stability of the acceptance in the Collins and Sivers angles 

A stringent test on the $ 3 (j = C, s) dependence of the acceptance ratio, already used for 
the previously published data, consists in checking the $j dependence of the ratios: 

R(t>) _ KM ■ KM ,261 

RA ^ ~ KM) ■ KM) m 

At the first order, assuming the absolute value of the target polarisation to be the same in each 
cell before and after reversal, it is 

In the very likely case in which af u ($j)/ajj($j) = <i f JA\j) : '<i - J. ^/) it is 



F + 
r d 



F + ■ F~ 




2 



(28) 



thus, the constancy in $j of Rj(Q>j) implies for each cell the ratio of the acceptances before 
and after the reversal to be constant in It must be noted that this constancy is not required 
to have unbiased estimators in the case the asymmetries are evaluated with the RPM. Still, this 
test is quite convincing on the stability of the apparatus. 

The ratios Rj were calculated in each bin of x, z, p\, for the Collins and Sivers angles, 
for leading and all hadrons and for both charges. In each bin the $ 3 distribution was fitted with 
a constant. Fig. [13] shows an example of the R c values for the Collins leading hadron sample 
in the 9 s-bins for the first period of data taking in 2004. The lines are the results of the fit. The 
quality of these fits are very good, as can be seen from Fig. fl4l where the distribution of the 
X 2 of all the fits is compared with the expected x 2 distribution for v — 15 degrees of freedom 
normalised to the number of entries. To conclude, this test gave an independent indication on 
the stability of the acceptances. 

Summarising, all the test we performed on the effects of acceptance and on fluctuation in 
time of the measured asymmetries gave results statistically compatible with what expected in 
the absence of systematic effects. 

5.4 Further tests on the fit of the ratio product quantities 

The quality of the fit of the double ratio quantities A™ with the function p ■ (1. + A™ ■ 
sin$j) has been checked looking at the x 2 distribution. Fig. [15] shows the x 2 distribution for 
Collins and Sivers asymmetries, in each bin of x, z, p\ for leading hadrons of both charges 
and for all the 5 data taking periods in 2002, 2003, and 2004. The curve is the expected x 2 
distribution for 14 degrees of freedom normalised to the number of entries. A perfect agreement 
can be observed. 

The effect of using two other different fitting functions, namely (1 + A m j ■ sin$j) and 
(po + A m j ■ sin as well as changing the $j bin size (8 bins instead of 16) have also been 
investigated. In both cases the effects turned out to be negligible, of the order of a few percents 
of the statistical error. 
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Figure 13: Distribution of the Rc values for the Collins leading hadron sample vs. x for the first 
period of data taking in 2004. The line shows the result the fit with a constant value. 



5.5 Different estimators 

To estimate the size of possible systematic effects, the asymmetries have also been eval- 
uated using two other estimators which rely on different assumptions of the acceptance varia- 
tions, namely the "difference method" (DM, in the following) and the "geometric mean method" 
(GMM), which are commonly used in evaluating asymmetries. 

The DM was used to extract the final values of the asymmetries published in [15]. The 
evaluation of the asymmetries is performed separately for the two target cells and, after checking 
their compatibility, they are combined by taking weighted averages. The asymmetries ec,k and 
es,k, where k = u,d indicates the target cell, were evaluated from the number of events with the 
two target spin orientations (+ spin up, and - spin down) by fitting the quantities 



J 



as 



(29) 



with the functions ec,k • sin $c and es,k ■ sin $5. The normalisation factor was taken equal to 
the ratio of the total number of detected events in the two orientations of the target polarisation: 
r = N^ tot k /N^ tot k . This procedure is correct if the difference of the absolute value of the target 
polarisation before and after the spin reversal is negligible. For what concerns the acceptances, 
they cancel in Eq. (l29l as long as the ratio a^ k (^j) /aJ k (Q>j) is constant in Qj. 

In the GMM, again the asymmetries are evaluated separately for the two target cells and 
then combined by taking weighted averages. The method consists in building the measured 
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quantities 

G&(*;) = V , — V y — - (30) 

Assuming: 

- a negligible difference in the absolute value of the target polarisation before and after the spin 
reversal (as for the DM), and 



CM 




Figure 15: x 2 distribution of the two parameter fit used to extract the asymmetries A™ compared 
to the normalised x 2 distribution for ndf = 14. 
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- that the acceptances satisfy the relation 

_ gZfcW (31) 

the asymmetries are evaluated by fitting the quantities G™ k ($j) with the functions ec,k • sin $c 
and £5^. • sin $5 in the range < $c,s < 7T- 

The advantage of this method is that the luminosity cancels. Still, the requirement on the 
acceptance stability is more demanding than in the case of the ratio product method, which was 
thus used to evaluate the final asymmetries. 

The fitted asymmetry values from the three methods were very close. Given the advan- 
tages of the RPM, the evaluation of the asymmetries with the GMM and DM has been consid- 
ered only a cross-check of the result, and the results of the experiment given in the following 
are those obtained using the RPM method. 



5.6 2-D fits 

In an entirely different approach, we have estimated the Collins and Sivers asymmetries 
using the standard linear least square method (LSM) and fitted in each kinematics bin the mea- 
sured asymmetries 

Aj (<f> h , <j) S ) = _ 77 7 v • JL 77 7 x (32) 

N j, u \<t>hAs) N j4 {(j) h ,<f)s) 

with the function 

H(<f> h , 0s) = ai + a 2 sin(0 /l + 5 - vr) + a 3 sin(30 ft - (f> s ) + 

+a 4 sm((f) h - <j> s ) + a 5 cos(0 ft - <j> 3 ) . (33) 

The first parameter should be one. The second, third and fourth terms arise from the 
transverse component of the target nucleon spin when the lepton beam is unpolarised, while 
the last term originates from the interaction between a longitudinally polarised lepton beam and 
a transversely polarised target (see ref [42]). The second term is the Collins term. The fourth 
term is the Sivers term. The last term has also physics interest of its own, it is related to the g 1T 
transverse momentum distribution function. 

In order to have enough statistics in each (<fi h , <p s ) bin, we have plotted the data in 8 bins 
of 4> h and 8 bins of 4>s- Each kinematic bin is thus split in 64 (<p h , <j> s ) bins. We have performed 
both 5-parameter fits and 4-parameter fits, in this second case fixing al to its expected value of 
1. Also, we have done two independent fits, one minimising a x 2 with a linear LSM, the second 
using MINUIT [43]. The results of the four fits are in excellent agreement, and both the fitted 
Collins and Sivers asymmetries and their errors turn out to be essentially identical to the values 
given by the one dimensional fits, as expected from the orthogonality of the different terms. The 
Collins and Sivers asymmetries as given by the 2-D fit turn out to be slightly correlated (the 
correlation coefficient ranges from -0.25 to 0.25), a known effect due to the considerably non 
uniform population of the 64 (0^, <f>s) bins. 

A full discussion of the procedure and of the results will be the subject of a separate paper. 



5.7 Systematic errors 

As spelled out in the previous sections, all tests performed in the different data taking 
periods did not give any evidence for the presence of systematic effects. 
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Figure 16: Distributions of Pkj (see text) for the leading hadron (left) and all hadron (right) 
asymmetries. 



The conclusive test was to look at the compatibility of the physics results obtained sep- 
arately for all the data taking periods. As already mentioned, from 2002 to 2004 data were 
collected in 5 periods and the compatibility test is significant. It has to be stressed that the 
2002 data from which the results have already been published, have been reanalysed using the 
slightly different event selection and analysis described here. The published data turned out to 
be very close to the new ones, with differences of the order of half of the systematic errors 
only in the less populated bins, essentially because of the different method used to evaluate the 
asymmetries (RPM instead of the DM used previously). 

The test was performed separately for all hadron and leading hadrons asymmetries, fol- 
lowing the procedure described in Section I5T21 The Collins and Sivers asymmetries in each x, 
Pt and z bin, for positive and negative hadrons, obtained in the 5 periods were compared with 
their weighted mean. The pulls are defined as 



Am 



A? 



a, 



1,5 



(34) 



where AJ 1 are the weighted means of the asymmetries, and are expected to follow a normal 
standard distribution. In total there are 500 Pkj values for the leading hadron asymmetries and 
520 for the all hadron asymmetries. 

The distribution of the Pkj values for all hadrons and for leading hadrons are shown in 
figure [161 As expected, these pulls follow a normal distribution. A very good agreement with 
purely statistical fluctuations has been seen for all the various subsamples of the asymmetries; 
the RMS of the different distributions are given for completeness in table |4j 

Since even in this last test we could not observe any indication for systematic effects, we 
concluded that the systematic errors due to acceptance and efficiency effects are considerably 
smaller than the statistical errors. 

The asymmetry scale uncertainty due to the uncertainties on Pt is of 5%. The error on 
the dilution factor /, which takes into account the uncertainty on the target composition, is of 
the order of 6%. When combined in quadrature, these errors give a global scale uncertainty of 
8%. 
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Table 4: RMS of the pull distributions for the different samples of asymmetries 



Asymmetries 


RMS 


positive hadrons 


leading 


Collins 
Sivers 
Collins 
Sivers 


0.882 ±0.056 

0.895 ±0.057 
agio -un n^7 

0.974 ±0.060 


all 


negative hadrons 


leading 


Collins 
Sivers 
Collins 
Sivers 


i rvn _i_n rififi 
0.971 ±0.061 
1.092 ±0.068 
0.991 ±0.061 


all 


leading 
all 


0.950 ±0.030 
0.996 ±0.031 


positive 
negative 


0.919 ±0.029 
1.026 ±0.032 



6 Results and comparison with models 

6.1 Measured asymmetries 

The results from the different data taking periods have been combined by making the 
standard weighted mean. 

Plots of the measured values of the Collins and Sivers asymmetries A Co ii and A Siv for the 
2003-2004 leading hadron sample against the three kinematic variables x, z and are given 
in Fig.[T7] Full points and open points refer to positive and negative hadrons respectively. The 
errors shown in the figure are only statistical. The same asymmetries are shown in Fig. [TS] for 
the 2003-2004 all hadron sample. Again, the errors shown in the figure are only statistical. 

The improvement in statistics with respect to the 2002 published result is clearly visible 
in Fig. [l9j where, as an example, the published Collins asymmetry data for all positive (left) 
and all negative (right) hadrons are compared with the 2003-2004 results. 

As already mentioned, the 2002 data have been reanalysed using the event selection and 
the analysis described in this paper. The effect on the measured asymmetries is very small, and 
is mainly due to use of the RPM in evaluating the raw asymmetries. As an example, Fig. |20| 
shows the comparison between the published and the new results for the Collins asymmetry. 
The new values from the 2002 data have been combined with the results from the 2003-2004 
data to evaluate the final asymmetries. 

The overall results from 2002-2004 deuteron target for the leading hadron sample and for 
the all hadron sample are given in Fig.|^T]andE2respectively. All these measured asymmetries 
are available on HEPDATA [41]. 

6.2 Comments and comparison with models 

As apparent from Fig. |2l]and|22j all the measured asymmetries are small, if any, and 
compatible with zero. This trend already characterised the published data of the 2002 run, and 
is confirmed by the new data with considerably improved precision. Small asymmetries is not 
a surprise. From the very beginning it was predicted that transverse spin effects be small in 
the deuteron due to the opposite sign which was expected for the u and d distributions, causing 
cancellations on the asymmetries of an isoscalar target, very much like in the helicity case. Still, 
it was not obvious that they would have been so small. 

An analysis of the results on the deuteron can be done only in conjunction with corre- 
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Figure 17: Collins asymmetry (top) and Sivers asymmetry (bottom) against x, z and f° r 
positive (full circles) and negative leading hadrons (open circles) from 2003-2004 data. Error 
bars are statistical only. In all the plots the open circles are slightly shifted horizontally with 
respect to the measured value. 



sponding proton data, which up to now have been measured only by the HERMES Collabo- 
ration [14]. The measured non zero Collins asymmetry on the proton has provided convincing 
evidence that both the transversity distribution A T u(x) and the Collins mechanism A° D^(z) 
are not zero. Independent evidence that the Collins mechanism is a real measurable effect has 
come from the recent analysis of the BELLE Collaboration [44,45]. Furthermore, the HERMES 
data on a proton target have provided convincing evidence that the Sivers mechanism is also at 
work. 

It is fair to say that the accuracy of the present HERMES data has allowed to extract only 
the leading contribution to the proton transverse asymmetry, that is the u quark contribution. 
Also, the interpretation of the HERMES data has led to surprising assumptions about the relative 
size of the favoured and unfavoured spin dependent fragmentation functions. A global analysis 
using all the available deuteron and proton data should allow now to provide first estimates of 
both the u and the d quark contribution, and clearly constitutes the next step in this work. 

In the following, first naive expectations, based on the simple formulas of Sections 11.21 
and ll.31 are given for the deuteron asymmetries, then the new deuteron data are compared to a 
few existing model calculations. 
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Figure 18: Collins asymmetry (top) and Sivers asymmetry (bottom) against x, z and f° r 
positive (full circles) and negative hadrons (open circles) from 2003-2004 data. Error bars are 
statistical only. In all the plots the open circles are slightly shifted horizontally with respect to 
the measured value. 



6.2.1 Collins asymmetry 

Although the measured deuteron asymmetries refer to unidentified hadrons, in the fol- 
lowing it will be assumed that the hadrons be pions (actually more than 80% are pions), so 
that the algebra considerably simplifies. Further simplification can be obtained by neglecting 
the sea contribution (i.e. A T q = A T s = and q = s = 0) and considering only the range 
0.1 < x < 0.3. This is justified by the fact that the PDFs are expected to be considerably differ- 
ent from zero in the valence region, and the HERMES data show non-zero values in the range 
0.05 < x < 0.3. Assuming 

Df = Dl = A , D? = Dl = D 2 , (35) 
A° Dl + = A° Df = A a T D x , A° T Df = A° D£ = A° T D 2 , (36) 

and using Eq. ©, one gets for n + on a proton target: 



+ AA T u v A^Di + A T d v A° T D 2 
CM ^ 4u v D 1 + d v D 2 (J/) 



and for tr : 



4A T u v A^D 2 + A T d v A^D 1 
Au v D 2 + d v D 1 



A p 6lu * ^^^ T ^^' . (38) 
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Figure 19: Collins asymmetry against x for all positive (left) and all negative (right) hadrons 
from 2002 data [15] (triangles) and from 2003-2004 data (circles). The 2002 values are slightly 
shifted horizontally with respect to the measured value. 

If, as suggested from the HERMES data, A° A = -A° D 2 , and taking D 2 ~ 0.5£>i, d v ~ 
0.5u v , the previous expressions become 



, A T u v a° d 1 

and 



AcoU ~ 2.5 u v D 1 (4U) 




Figure 20: Collins asymmetry against x for positive (left) and negative (right) hadrons from 
2002 data. The triangles are the published results [15] and the circles the results of the new 
analysis. The published values are slightly shifted horizontally with respect to the measured 
value. 
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Figure 21: Overall results for Collins asymmetry (top) and Sivers asymmetry (bottom) against 
x, z and p\ for positive (full circles) and negative leading hadrons (open circles) from 2002, 
2003, and 2004 data. Error bars are statistical only. In all the plots the open circles are slightly 
shifted horizontally with respect to the measured value. 



respectively. As already stressed, the u-quark contribution is dominant. 

For a deuteron target equations (l37l) . (l38l) . (l39t . and (l40l) become respectively 



A 



A T u v + A T d v AA° T D X + A° T D 2 



Coll 



d, , ■ 



a 



(,7T 

Coll 



A T u v + A T d v A%D X + AAtD 2 



D 1 + 4L> 2 



and 



A ~ 
^Coll — 



^Coll — 



3 A T u v + A T d v A^Dx 



7 u v D% 
3 A T u v + A T d v A° T Dx 



4.5 



u. 



(41) 



(42) 



(43) 



(44) 



Both the 7r + and the n~ Collins asymmetries on the deuteron are proportional to A T u v (x) + 
A T d v (x), therefore cancellation is expected to reduce considerably the effect which has been 
measured on the proton. As a matter of fact, assuming A T d v = (no limit on the size of A T d v is 
provided by the HERMES data) one derives relations between the Collins asymmetry measured 
by HERMES and COMPASS which are only marginally satisfied. Thus, the present precise data 
on Aqou should allow to extract Ard v . 
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Figure 22: Overall results for Collins asymmetry (top) and Sivers asymmetry (bottom) against 
x, z and p\ for all positive (full circles) and all negative hadrons (open circles) from 2002, 2003, 
and 2004 data. Error bars are statistical only. In all the plots the open circles are slightly shifted 
horizontally with respect to the measured value. 



This was not the case in so far. Three global analyses have been performed with the 
published data, trying to derive bounds on the transversity distributions and the Collins frag- 
mentation functions. In Ref. [46] the Soffer bound |Ayg| — (q + Ag)/2 was used, a fit of the 
HERMES data set was performed, and the Collins functions were extracted. Two different sce- 
narios for favoured and unfavoured Collins fragmentation functions were considered, but the 
fits always favoured a relation A^D 1 ~ — A° D 2 . The comparison of the fit results with the 
COMPASS data shows a fair agreement, as apparent from Fig. |23l although the data do not 
exhibit the trend with x which is suggested by the model. The upper and lower curves in the 
figures correspond to the 1-cr errors of the fitted parameters. 

In Ref. [27] a chiral quark-soliton model was used for the transversity distributions, and 
the Collins fragmentation function was derived from a fit to the HERMES data, which do not 
constrain the Axd distribution. A comparison with the present COMPASS results shows again 
a fair agreement (Fig.l2"4l. The upper and lower curves in the figures correspond to the uncer- 
tainty in the Collins fragmentation functions as obtained from the fit. Independent extraction 
of the Collins function was performed by fitting the BELLE data. The result was found to be 
compatible with the one obtained fitting the HERMES data. 

Similar results were obtained in Ref. [47]. Two different scenarios were used for transver- 
sity, either A^g = Ag, or the Soffer bound, and the Collins fragmentation functions were ex- 
tracted from a fit to the HERMES data. The fits were very good in both cases. The extracted 
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Figure 23: Comparison between the present results for positive (left) and negative (right) 
hadrons with the calculations of Ref. [46] (scenario 1). 



Collins functions were then used to reproduce the published 2002 COMPASS data. The agree- 
ment is acceptable for both scenarios, as apparent from Fig. |25J although also in this case the 
expected increase of A Co u with increasing x is not manifested by the COMPASS data. The up- 
per and lower curves in all the figures correspond to 1-er deviations of the parameters. Also in 
this case the BELLE data are reproduced, and one can conclude that the Collins mechanism in 
SIDIS and in e + e~ — > hadrons are the same. 

To summarise, the new data are compatible with the uncertainty bands given by the phe- 
nomenological calculations, but the trend of the data with x is not the one suggested by the 
central value of the calculations: this is a clear indication that the COMPASS data will surely 
help in constraining the parameters of the models. 
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Figure 24: Comparison between the present results for positive (left) and negative (right) 
hadrons with the calculations of Ref. [27]. 
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Figure 25: Comparison between the present results for positive (left) and negative (right) 
hadrons with the calculations of Ref. [47], where two scenarios were tested for transversity: 
A T q = Aq was assumed in one scenario (upper figures), while A T q — (q + Aq)/2 was as- 
sumed for the calculations of the lower figures. 

6.2.2 Sivers asymmetry 

Also in this case it is useful to consider the expressions one obtains for A Siv in the hy- 
pothesis that all hadrons are pions. Again, the simplified analysis is restricted to the valence 
region. 

Neglecting the sea contribution (i.e. Aq q = Aq s = and q = s = at all x) and 
assuming Dl + = = D\ and D% + = = D 2 , on a proton target, from Eq. (fl4l one gets 
for 7r + : 

, i7r + _ . 4A^£> 1 + A%d v D 2 



A Siv 



4u v D 1 + d v D 2 



(45) 



and for tt 



- _ AA^u v D 2 + A^d v D t 
Siv ~ \u v D 2 + d v D x 



(46) 



Assuming D 2 ~ 0.5Di, d v ~ 0.5u v , the previous expressions become 



^Siv — 



A%u v 



(47) 
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and 

2A T v + A r d 



respectively. Since the Sivers asymmetries for n as measured by HERMES is about zero, in 
this very simplified treatment it follows that 

Ajfo = -2A£u„ . (49) 

For a deuteron target the Sivers asymmetries can be written as 

& - (50) 



and 



which implies A d ^ ~ Ag^ . The approximatively zero Sivers asymmetries for positive and 
negative hadrons observed in COMPASS require 

A T Q d v ~ -A? Uv , (52) 

a relation which is also obtained in some models. 

From Eq. (l47l) . d48t . and (1491 . a relation between the Sivers asymmetry measured by 
COMPASS and those measured by HERMES can be derived 

A P,7T+ 

which is only marginally satisfied by the data. 

Also the Sivers data have been looked upon independently by three different groups. 
In Ref. [26] a fit of the HERMES data was performed, on the assumption that A^d(x) = 
—Aqu(x). A good agreement with the HERMES data was obtained, and a zero asymmetry in 
case of a deuterium target was predicted, as shown in Fig. [26] The curves in the figure indicate 
the expected size of the effect on A$i V of the l/iV c -corrections. The COMPASS data fall well 
within the band resulting from the model. 

The authors of Ref. [48] could estimate the Sivers functions by fitting both the HERMES 
and the published COMPASS data getting A%d v ~ -A%u v . Leading Order MRST01 sets of 
unpolarised distribution functions [49] were used, together with Kretzer's set of fragmentation 
functions [50]. A very good agreement with the experimental data was reached, as apparent 
from Fig.|27j where the upper and lower curves correspond to l-a deviation of the parameters. 
In this model the z and "p\ dependence of the single spin asymmetries are also well described. 

In Ref. [46] it was assumed that the final hadron transverse momentum is the transverse 
momentum in the Sivers function, i.e. a collinear fragmentation was assumed. GRV98 leading 
order distribution functions [51] were used, along with Kretzer's fragmentation functions. The 
fit of the HERMES data is very good, and gave as a result A$d v ~ — Aq u v , but the prediction 
for COMPASS are not in agreement with the new data, as shown in Fig. The upper and 
lower curves in the figures correspond to l-a errors in the fitted parameters. Here again, a new 
global fit using the deuteron data should be able to extract the Sivers function for the d-quark 
and improve the agreement with the data. 

Very recently, the smallness of the Sivers asymmetry for positive and negative hadrons 
on the deuteron has been interpreted as evidence for the absence of gluon orbital angular mo- 
mentum in the nucleon [52]. 
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Figure 26: Comparison between the present results for positive (left) and negative (right) 
hadrons with the calculations of Ref. [26]. 



7 Conclusions 

After providing the very first SIDIS data on transverse spin asymmetries on a transversely 
polarised deuteron target, COMPASS now publishes the overall results from the data collected 
in 2002, 2003, and 2004, increasing the statistics as compared to the published 2002 data by 
a factor of ~7, so that even at large x the errors on the measured asymmetries are only a few 
percent. 

All the measured asymmetries are small, mostly compatible with zero within the mea- 
surement errors. Presently, the most likely interpretation, taking into account the correspond- 
ing measurements of the HERMES collaboration on a proton target, is that in the COMPASS 
isoscalar target there is a cancellation between the proton and the neutron asymmetries. Also, 




Figure 27: Comparison between the present results for positive (left) and negative (right) 
hadrons with the calculations of Ref. [48]. 
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Figure 28: Comparison between the present results for positive (left) and negative (right) 
hadrons with the calculations of Ref . [46] . 

the independent evidence provided by the BELLE data that the Collins effect is a real physi- 
cal mechanism guarantees that the transversity distributions A T q(x) can be extracted from the 
single spin asymmetries measured in SIDIS. A global analysis of all the presently available 
data is now mandatory, and the inclusion of the new precise deuteron data from COMPASS will 
surely allow to provide first estimates of the down quark transversity distribution A T d and of the 
Collins fragmentation functions A^D. Already now the COMPASS data for the Sivers asym- 
metry provide convincing evidence on the cancellation of the u and d quark Sivers distribution 
functions. 

It has to be stressed, however, that the measured effects are rather small, of the order of 
a few per cent at most, and that flavour separation requires much larger statistics than presently 
available. Also, within the present errors, the physics interpretation is not straightforward and 
has led to some surprises. The present data are of fundamental importance because they have 
opened up the road to transverse momentum dependent distribution and fragmentation func- 
tions, but they will by no means suffice to determine these new functions. New data are needed, 
and particularly, new proton data. HERMES is still finalising the analysis of their 2005 proton 
run, which will double the statistics of their analysed data. The BELLE Collaboration is pro- 
ducing more and more accurate data on the Collins fragmentation function. In the near future, 
COMPASS plans proton runs, which should result in particularly reduced error bars at large x, 
thanks to the much improved geometrical acceptance of the spectrometer which is obtained us- 
ing the new COMPASS polarised target magnet. A global analysis will then again be necessary, 
and from all those measurements it should be able to provide the u and d quark transversity 
distributions. 



8 Acknowledgement 

We thank M. Anselmino, X. Artru, V. Barone, and A. Prokudin for many interesting 
discussions. We gratefully acknowledge the support of the CERN management and staff, as 
well as the skills and efforts of the technicians of the collaborating institutes. 



38 



References 

[1] G. Bunce etal, Phys. Rev. Lett. 36 (1976) 1113. 
[2] G.L. Kane etal., Phys. Rev. Lett. 41 (1978) 1689. 

[3] see, f.i., A. Efremov and O. Teryaev, Sov. J. Nucl. Phys. 39 (1984) 962, and the papers of 

the Lund team on A polarisation. 
[4] J. Antille etal., Phys. Lett. B 94 (1980) 523. 

[5] V.D. Apokin et al, in 'Proceedings of the 7-th International Symposium on High En- 
ergy Spin Physics, Protvino, USSR, 1986 (Institute of High Energy Physics, Serpukhov, 
USSR, 1987) Vol 2, pag. 83. 
[6] The Fermilab E704 Collaboration, A. Bravar et al, Phys. Rev. Lett. 77(1996)2626. 
[7] European Muon Collaboration, J. Ashman et al, Phys. Lett. B 206 (1988) 364; Nucl. 

Phys. B 328 (1989) 1. 
[8] J.P Ralston and D.E. Soper, Nucl. Phys. B 152 (1979) 109. 
[9] X. Artru and M. Mekhfi, Z. Phys. C 45 (1990) 669. 
[10] R.L. Jaffe and X. Ji, Phys. Rev. Lett. 67(1991) 552. 
[11] J. Collins, Nucl. Phys. B 396 (1993) 161. 
[12] G. Bunce etal, Part. World 3 (1993) 1. 

[13] The HELP Collaboration, B. Vuaridel et al, CERN reports LEPC 93-14, LEPC/P7 
(1993). 

[14] The HERMES Collaboration, A. Airapetian et al, Phys. Rev. Lett. 94 (2005) 012002. 

[15] The COMPASS Collaboration, V. Y. Alexakhin et al, Phys. Rev. Lett. 94 (2005) 202002. 

[16] V. Barone and PG. Ratcliffe, Phys. Rept. 359(2002) 1. 

[17] R.L. Jaffe and X. Ji, Nucl. Phys. B 375 (1992) 527. 

[18] A. Kotzinian, Nucl. Phys. B 441 (1995) 234. 

[19] P.J. Mulders and R.D. Tangerman, Nucl. Phys. B 461 (1996) 197. 

[20] X. Artru, "The Transverse Spin", Proceedings of 10th Rhodanien Seminar, 

|hep-ph/0207309l 
[21] D. Sivers, Phys. Rev. D 41 (1990) 83 

[22] S. J. Brodsky, D. S. Hwang and I. Schmidt, Phys. Lett. B 553 (2003) 223. 

[23] J. C. Collins, Phys. Lett. B 536 (2002) 43. 

[24] A. V. Belitsky, X. Ji and F. Yuan, Nucl. Phys. B 656 (2003) 165. 

[25] D. Boer and P.J. Mulders, Phys. Rev. D 57 (1998) 5780. 

[26] J. C. Collins et al, Phys. Rev. D 73 (2006) 014021. 

[27] A. Efremov et al, Phys. Rev. D73 (2006) 094025. 

[28] The COMPASS Collaboration, E. S. Ageev et al, Phys. Lett. B 633 (2006) 25. 

[29] The COMPASS Collaboration, E. S. Ageev et al, Phys. Lett. B 612 (2005) 154; the 
COMPASS Collaboration, V. Y. Alexakhin et al, "The Deuteron Spin Dependent Struc- 
ture Function gf and its First Moment", in preparation. 

[30] The COMPASS Collaboration, "The COMPASS experiment at CERN", in preparation. 

[31] J. Ball et al, Nucl. Inst. Meth. A 498 (2003) 101. 

[32] The Spin Muon Collaboration, D. Adams et al, Nucl. Instrum. Meth. A 437 (1999) 23. 
[33] Y. Giomataris et al, Nucl. Inst. Meth. A 376 (1996) 29; F. Kunne, Nucl. Phys. A 721 
(2003) 1087c. 

[34] F. Sauli, Nucl. Inst. Meth. A 386 (1997) 531; B. Ketzer etal, IEEE Trans. Nucl. Sci. 49 
(2002) 2403. 

[35] V.N. Bychkov et al, Nucl. Inst, and Meth. A 556 (2006) 66. 
[36] E. Albrecht et al, Nucl. Inst. Meth. A 504 (354) 2003. 



39 



[37] C. Bernet et al, Nucl. Inst. Meth. A 550 (2005) 217. 
[38] L. Schmitt etal, IEEE Trans. Nucl. Sci. 51 (439) 2004. 
[39] for all the COMPASS off-line pakages see 

http://wwwcompass.cern.ch/compass/software/offline/welcome. html#software. 
[40] R. Brun and F. Rademakers, Nucl. Instrum. Meth. A 389 (1997) 81. See also 

http : //root . cern . ch/. 
[41] HEPDATA: The Durham HEP Database, 

http://durpdg.dur.ac.uk/HEPDATA/index.html. 
[42] A.V. Efremov, Annalen Phys. 13 (2004) 651. 
[43] F. James, MINUIT, CERN Program Library Long Writeup D506, 

http://wwwasdoc.web.cern.ch/wwwasdoc/WWW/minuit/minmain/minmain.html. 
[44] The Belle Collaboration, K. Abe et al, Phys. Rev. Lett. 96 (2006) 232002. 
[45] A. Ogawa et al, "Spin dependent fragmentation function at Belle," 

|arXiv:hep-ex/06070T4| 
[46] W. Vogelsang and FYuan, Phys. Rev. D72 (2005) 054028. 
[47] M. Anselmino et al, RICH + AGS annual Users' Meeting June 6, 2006, 

www.bnl.gov/rhic_ags/users_meetings/workshops/basp, and paper in preparation. 
[48] M. Anselmino et al, Phys. Rev. D 72 (2005) 094007. 
[49] A. D. Martin et al, Phys. Lett. B 531 (2002) 216. 

[50] S. Kretzer, "Fragmentation functions from flavour-inclusive and flavour-tagged e+ e- 

Phys. Rev. D 62 (2000) 054001. 
[51] M. Gliick, E. Reya and A. Vogt, Z. Phys. C 67 (1995) 433. 

[52] S. J. Brodsky and S. Gardner, "Evidence for the Absence of Gluon Orbital Angular Mo- 
mentum in the Nucleon", arXiv:hep-ph/0608219. 



40 



